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Abstract 

In this paper we address the adequacy of various approximate methods of 
including Coulomb distortion effects in (e, e') reactions by comparing to an 
exact treatment using Dirac-Coulomb distorted waves. In particular, we ex¬ 
amine approximate methods and analyses of (e, e') reactions developed by 
Train! et al. using a high energy approximation of the distorted waves and 
phase shifts due to Lenz and Rosenfelder. This approximation has been used 
in the separation of longitudinal and transverse structure functions in a num¬ 
ber of (e, eb experiments including the newly published ‘^^^Pb{e, e') data from 
Saclay. We find that the assumptions used by Traini and others are not valid 
for typical (e, e') experiments on medium and heavy nuclei, and hence the 
extracted structure functions based on this formalism are not reliable. We 
describe an improved approximation which is also based on the high energy 
approximation of Lenz and Rosenfelder and the analyses of Knoll and com¬ 
pare our results to the Saclay data. At each step of our analyses we compare 
our approximate results to the exact distorted wave results and can there¬ 
fore quantify the errors made by our approximations. We hnd that for light 
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nuclei, we can get an excellent treatment of Coulomb distortion effects on 
(e, e') reactions just by using a good approximation to the distorted waves, 
but for medium and heavy nuclei simple additional ad — hoc factors needs 
to be included. We describe an explicit procedure for using our approximate 
analyses to extract so-called longitudinal and transverse structure functions 

from (e, e') reactions in the quasielastic region. 
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I. INTRODUCTION 


Medium and high energy electron scattering has long been acknowledged as a useful tool 
in the investigation of nuclear structure and nuclear properties, especially in the quasielastic 
region. In the plane wave Born approximation (PWBA), where electrons are described as 
Dirac plane waves, the cross section for inclusive quasielastic (e, e') processes can be written 
simply as 


^ u.) + |tau^ I - ^]Sr{q ,»)} 


( 1 ) 


where q't = is the four-momentum transfer, (Tm is the Mott cross section given by 


C.2 0 


(Tm — 


' 2E ' sm‘* ^ 


f, and Sl and St are the longitudinal and transverse structure functions 


which depend only on the momentnm transfer q and the energy transfer uj. By keeping the 
momentnm and energy transfers hxed while varying the electron energy E and scattering 
angle de-, it is possible to extract the two structure functions with two measurements. How¬ 
ever, when the electron wavefunctions are not Dirac plane waves, but rather are distorted by 
the static Coulomb field of the target nucleus, such a simple formulation as given in Eq. (1) 
is no longer possible and, in general, the cross section does not separate into the snm of 
longitndinal and transverse structure functions with coefficients which only depend on the 
electron kinematics. Clearly the size of the Coulomb distortion effects depends on the charge 
of the target nncleus and on the energy of the electrons. In general. Coulomb effects are not 
too large near the peaks of cross sections, but can have greatly magnified effects when one 
extracts “structure functions” by subtracting one cross section from another. For example, 
in a recent distorted wave calculation [1] of ^®0(e,e'p) in the qnasielastic region, we found 
Coulomb effects on the extracted spectroscopic factors to be approximately 3%, while the 
effects on the extracted fourth structure function was approximately 15%. 

However, some sort of extraction of strnctnre functions, albeit approximate, is still very 
appealing since structure functions are sensitive to different aspects of the underlying knock¬ 
out process or the hnal state interaction. Furthermore, (e, e') reactions in the qnasielastic 
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region are particularly appealing in that the cross section is the sum of many knockout pro¬ 
cesses and some of the various Coulomb effects may be partially averaged out. In fact, this is 
the case. A rather simple approximation known as the Effective Momentum Approximation 
(EMA) where the electron momenta pij are modihed by the value of the Coulomb potential 
at the center of the nucleus goes quite far in reproducing the Coulomb distortion effects for 
light nuclei. However, the EMA for heavy nuclei does not adequately reproduce the DWBA 
cross section in the quasielastic region. 

During the past decade, (e, e') cross sections have been measured [2]- [4] for a number 
of nuclei in the quasielastic region and either plane wave or EMA was used to extract 
longitudinal and transverse structure functions. These extracted structure functions were 
compared to the predictions of a simple Fermi-gas model, and in some cases, there appeared 
to be large suppression (up to about 40%) of the longitudinal structure functions. There were 
also disagreements with the extracted transverse structure functions and the predictions of 
the Fermi-gas model, but these were expected since exchange currents, pion production and 
other processes that are primarily induced by transverse photons were not included in the 
Fermi-gas model. It should be noted that the Fermi-gas model is a rather crude description 
of a nucleus. In particular, the shape of the structure functions at hxed momentum transfer 
plotted as a function of energy transfer is not well described. On the other hand, we found 
that a “single-particle” relativistic model using relativistic Hartree wavefunctions coupled 
with the full DWBA treatment of Coulomb distortion for the electrons was in good agreement 
with the measured cross sections in the quasielastic region for [5] and in reasonably 

good agreement for Furthermore, the longitudinal structure function extracted for 

was well produced in magnitude and shape by this model [6]. 

More recently, Jourdan [7] has examined the world data set for inclusive quasielastic 
scattering on and at momentum transfers of 300, 380 and 570 MeV/c, 

and has evaluated the Coulomb sum rule. No evidence of suppression is found for the 
highest q value (570 MeV/c) where the sum rule is most model independent. Jourdan 
used our Coulomb corrections in arriving at this conclusion. Thus, on the basis of our direct 
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comparison with the measured cross sections of and of Jourdan’s results, the published 
analyses [8] for e') which claim up to 50% suppression of the longitudinal structure 

function is surprising. However, the approximate treatment of Coulomb distortion used in 
the analysis is not accurate, and leads to doubts about the extracted structure functions. As 
a separate matter, we question the nuclear model used in making the claim of suppression. 

In this paper we investigate the possibility of including Coulomb distortion effects in 
(e, e') reactions in the quasielastic region for medium and heavy nuclei in an approximate 
way. We have an advantage as compared to previous workers in that we have an exact 
treatment of the static Coulomb distortion of the target nucleus via a Distorted Wave Born 
Approximation (DWBA) calculation to which we can compare [9]. In Section II we will 
discuss various approximations that permit a “plane-wave-like” approach to the treatment 
of Coulomb distortion and compare the approximate results to the exact DWBA results in 
a step by step way. We obtain an approximate potential due to the electron current which 
describes Coulomb distortion quite well. In Section III, we apply this potential with further 
approximations to the particular case of inclusive quasielastic processes. Finally we compare 
our calculations to the Saclay data [8] for the inclusive reaction e'). 

In addition, we give an explicit procedure for extracting longitudinal and transverse 
structure functions from inclusive cross section data in the quasielastic region from medium 
and heavy nuclei, and make some general conclusions. 

II. APPROXIMATIONS 

A. Approximation of the Electron Potential 

The four potential arising from the electron current is simply given in terms of the 
retarded Green’s function by 

A^(r) = J G{re,r)jf,{re)dre (2) 

where 
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G'(re,r) = 


giij|re-r| 


Tp — r 


and the electron cnrrent is given in terms of the initial and hnal electron wavefnnctions by 




( 3 ) 


For scattering processes this current extends over all space and thus the integral in Eq. (2) 
is not straightforward. 

Knoll proposed [10] a way of replacing the integral in Eq. (2) by a series of differential 
operators by using the transformation, 


I V{r - r')f{r')dr' = + *V)e-*'i''V(r) (4) 

between the function K(r — r') and its Fourier transform K(q'). Applying this result to 
Eq. (2) and making a Taylor series expansion, we obtain the following series expansion for 
the potential. 


= 3^ 




/2 I y2 




( 5 ) 


q'^ — q'^ — q'^ — 

Note that while the momentum variable q' is arbitrary, the choice affects the convergence 
of the series. In particular, for the case of Dirac plane waves for the electron the only 
dependence of the electron current on r is simply and choosing q' = q results in the 
vanishing of all the terms except the hrst, and Eq. (5) reduces to the well known Mdller 
potential. 


A^(r) = 

q^ — 

where u is the familiar Dirac plane wave spinor. 

As a test of this approximate procedure for calculating the potential, we calculated the 
electron charge distribution in the presence of the static Coulomb potential arising from 
the ground state charge distribution of using the partial wave solutions of the Dirac 

equation and evaluated the zeroth, zeroth plus hrst and zeroth plus hrst and second terms 
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in Eq. (5). Using these potentials we calculated the inelastic scattering cross section induced 
by a surface nuclear charge transition density 

Pn (r) = -^S{r - Rn)Y^{f) (6) 

tin 

where is the nuclear radius. The ground state density of the nucleus was described by 
a Fermi distribution of radius R = 6.65 fm and total charge Z = 82. The cross sections 
were calculated at initial energy Ei = 400 MeV and hnal energy Ef = 300 MeV with energy 
transfer uj = 100 MeV. In agreement with Knoll [10], we found that the first and second 
correction terms are sufficient to £11 up the minima and the contribution of the second 
correction term was less than 2% for momentum transfer q> 350 MeV/c. We conclude that 
this high momentum approximation provides an alternative procedure for calculating the 
potential arising from inelastic electron scattering processes. However, this procedure does 
require the numerical solution of the Dirac equation using partial waves which does take 
some computational time. One advantage, however, is that the radial functions only have 
to be calculated out to about three times the nuclear radius. 

We also applied this procedure to various approximate solutions of the Dirac equation and 
noted that the first and second order terms were not well controlled if one approximates the 
electron current. Clearly one should not approximate a function and then differentiate it. In 
the following when we examine approximate electron wavefunctions, and hence approximate 
electron currents, we will only use the zeroth term in Eq. (5), but will be sensitive to the 
choice of q'. 


B. High Energy Wavefunction Approximation 


The incoming Coulomb distorted electron scattering wave function that satisfies appro¬ 
priate boundary conditions for an electron with spin Si can be expressed in the form of a 
partial wave sum by [11] 




Si I 


= Ani 


lEi 


m 


2Ei 


E' 

K.,fl 




ll — Si Si 


yrfl — Si* 




( 7 ) 
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where the spinor ip^ is an eigenstate with angular momentum quantum numbers k and p 
given explicitly by 


where the spin-angle functions are 


igK{r)x-^{T) 


( 8 ) 




r = 


I 1/2 j 

/ u —5 5 /i 


Yf 


(f)x. 


( 9 ) 


The Dirac quantum number k determines the angular momentum labels for both I and 
j. Note that if we ignore the mass of the electron (valid for scattering angles away from 
extreme forward and backward angles), the electron helicity is a good quantum number, and 
we only require positive k solutions. The radial functions /«, and and their corresponding 
phase factors 6^ are obtained by numerically solving the Dirac radial equation for a hnite 
spherically symmetric nuclear charge distribution. The outgoing distorted wave function 
^lJy is found from Eq. (7) by making the replacements > —iSk), {si ^ Sf), {Ei Ef), 

and (pi —^ pj). Note that by setting the phases to zero and replacing /«, and g,^ by the 
spherical Bessel functions ji{pr) and sign{K,)ji{pr), where I = /(—«), the partial wave sum 
in Eq. (8) can be summed to give the Dirac plane wave solution. 


e*P", (10) 

<T p 

J 

where the electron spin label has been suppressed. Thus, one way of obtaining a plane-wave- 
like Coulomb distorted wavefunction, is to approximate the radial functions /«, and g^ by 
spherical Bessel functions, and to approximate the scattering phase shifts by an operator that 
can be pulled out of the partial wave sum. Using these two ideas, Lenz and Rosenfelder [12] 
obtained an approximate scattering wavefunction for high energy electrons which includes 
Coulomb distortion for a hnite nucleus in an approximate way more than twenty years ago. 

In the high energy limit with good helicity {E 3> me), the approximate wavefunction of 
Lenz and Rosenfelder [12], and Knoll [10] can be written as 


'hi(r) = 


IE, 


m 


2Ei 



( 11 ) 




where Up is the Dirac plane wave spinor and 

1 

rjir) = — \p — V{r')]dr' + correction term. (12) 

pr Jo 

Lenz and Rosenfelder approximated the dependence of the phases on the operator by 
(5(J^) = (5 i/2 + b(J‘^ — 3/4) where b depends on the Coulomb potential. This particular 
approximation for the phase shifts is not necessary since any function of the operator 
still allows the partial wave sum to be carried out. Lenz and Rosenfelder claimed that 
this equation is valid for j _|_ i/2<^pR. The second condition being primarily 

determined by the expression used for the phase shifts. 

The approximate radial functions are given by spherical Bessel functions with modihed 
argument, 


where x = p'{r)r and 


fnir) = —ji{x) 
pr 

X 

g^ir) = —sign{K)j-i{x) 
pr 


, 1 R A 

p{r)=p -/ V{r)dr-\ - 

r JO r 


(13) 


(14) 


where V(r) is the spherically symmetric Coulomb potential of the target nucleus and A is 

2 

a correction term in the argument of order and is normally neglected. 

To avoid having an r-dependent momentum, Traini et al. [14] made the further approx¬ 
imation that p'{r)=p'{0) = p — V{0) where R(0) is the static Coulomb potential evaluated 
at the origin. This approximation of the radial wavefunction leads to what is known as the 
Effective Momentum Approximation(EMA). An approximation to the radial wavefunction 
where we keep the r-dependence in the momentum but still neglect A will be referred to 
as the Local Effective Momentum Approximation(LEMA). In looking at these approximate 
wavefunctions, one should keep in mind that it is the spatial region around the nuclear 
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surface that contributes most significantly to electron induced transitions. That is, partial 
waves with angular momenta of order pR where R is the nuclear radius play a large role in the 
transition amplitnde. Thus, the approximation that the angular momentum is signihcantly 
less than pr is not valid. With this point in mind, we sought an ad — hoc approximation for 
the correction A. We found that the following expression describes the radial wavefunction 
at larger radial distances qnite well, 

1 fvfcf 

A{aZ, E, K^) = —ciZ{——Y = cbK^- 

tj 

The parameter a = —aZ{^Y where the nnmber 16 is given in MeV and was determined 
by comparing to the exact distorted radial wavefnnction. We will label this case LEMA + 
A. Note that we have chosen the parametrization of A so that it contains k? which can be 
expressed in terms of the operator so that we will still be able to sum the partial wave 
series. 

To investigate these approximations we compare the radial wavefnnctions calculated 
using EMA, LEMA, and LEMA + A to the exact Coulomb distorted waves (DW) for 
various angular momentum states for electron scattering from The radius of ‘^^^Rb is 

approximately 6.5 fm. Figs. (1) and (2) show the comparison between three approximate 
radial wavefnnctions and the exact distorted radial wavefunction for different energies {E = 
200 MeV, 400 MeV and 600 MeV) and different k values(fi; = 5 and k = 15). The EMA 
wavefnnction is acceptable at small radial distances, particularly much less than the nnclear 
radius, but at larger radii the approximate radial wavefunction is shifted too mnch to the 
left indicating the potential V (0) at the origin is too large. The LEMA wavefunction is 
a much better approximation to the exact wavefnnction than EMA for medium and high 
energies ont to radii beyond the nnclear radius, but at lower energies ( less than 250 MeV), 
it also deviates signihcantly from the DW resnlt. The wavefnnction with the LEMA +A 
correction agrees with the distorted one almost perfectly above 400 MeV. 
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FIG. 1. Radial wavefunctions in for k = 5. The top figure is for energy E = 200 MeV, 

the middle for energy E = 400 MeV and the bottom for energy E = 600 MeV. The solid line is 
the exact DW wavefunction while the dash-dotted line is LEMA-I-A, the dashed line is LEMA and 
the dotted line is EMA. 

At the lower energy this wavefunction is a bit to the right of the distorted wavefunction, 
but the discrepancy is acceptably small. Our conclusion is that LEMA is much better than 
EMA and may be acceptable for some reactions. The approximation LEMA -|- A furnishes 
quite a good description of the distorted wave radial functions out to several nuclear radii. 
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FIG. 2. Radial wavefunctios in 208 £qj. ^ same as Fig. (1). 

The other important ingredient in Coulomb distortion considerations are the phase shifts 
Sf^. In order to sum the partial wave series in Eq. (7), the phase shifts in Eq. (11) were 
expressed as a function of the total angular momentum operator J^. The eigenvalues for 
this operator are j{j + 1) when operating on the partial waves, so the question is how well 
are the exact phase shifts reproduced by the expression [10], 

Sj = Si/2 + b[j{j + 1) - -] 

= <^1/2 + — 1] (15) 
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where we used the relation j =\ k \ —1/2. For a uniform charge distribution of radius R, 





In 2pR) + b 


( 16 ) 


and 


3Za 

4p(0)'^R2- 


(17) 


We investigated this approximation for the phase shifts for the Coulomb potential with 
‘^^^Pb and 400 MeV electrons. The phase shifts are in good agreement with the exact phase 
shifts for small k values, but for large k values, the approximate phase shifts are much too 
large in magnitude. The breakdown of the approximation occurs for K^pP as expected. 
However, as noted earlier, it is these orbitals that play a dominant role in electron induced 
reactions from the nucleus. In order to avoid this violation, we assume that the approximate 
phase shifts contain a expressed as 


5^ — Oq T (I 2 K, T Qj^K 

= ho + b2[j{j + 1 ) - ^] + bi[j{j + 1 ) - (18) 

where 60 = Oo + 02 + ^ 4 , &2 = <^2 + 204 , and 64 = 04 . Unfortunately, we do not have a simple 
analytical expression for the coefficients, so Oq, 02 and 04 are htted to the exact phase shifts 
calculated with a distorted wave code. As is evident in Fig. (3), including the terms leads 
to a much better description of the exact phases. 
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FIG. 3. Comparison of the exact and fits to the phases in for the energy E = 400 

MeV and Kmax = 35. The diamonds are the exact phases and the dashed line is the K^-fit to the 
phases and the solid line includes the term in the fit. 

Clearly higher powers of could also be included if needed. Of course, using the 
description of the phases given in Eq. (18) rather than the one in Eq. (15) requires calculation 
of the exact Coulomb phases for both the incoming and outgoing electron energy and then 
determining the coefficients by a htting procedure. However, the solution for the phases 
is straightforward and very rapid for modern computers, so this poses no real practical 
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problem. 

We have calculated various multipoles of the scalar potential in the partial wave formal¬ 
ism and conhrmed that the approximate radial function and the approximate phases are in 
close agreement with the exact potential. These two approximations permit the summa¬ 
tion of the partial wave series of Eq. (7) so that the electron wavefunction for incoming or 
outgoing waves can be written as 


^(±)(r) 


f^g±i5(Jbghi+7^)p'h)-r^ 

p 


(19) 


The operator is the square of the total angular momentum operator (J = L + S), and 
previous workers [14,15] have made a power series expansion of the phase term and applied 
successive terms to the plane-wave-like part of the wavefunction. The problem with this 
procedure is that the phase shifts are greater than one for k values of importance and 
this series is very slowly converging. It is straightforward to check that keeping only the 
hrst three terms in the exponential expansion leads to signihcant errors. We choose not to 
make this expansion, but to approximate by the orbital angular momentum squared 
and further to replace by its classical value (rxp(r)')^. Clearly we are neglecting the 
spin dependence of the phase shifts by this classical approximation, but since the processes 
we are interested in are dominated by angular momentum values around 10 or more, we 
expect the spin dependence to be negligible. We conhrmed this estimate by comparisons 
of our partial wave calculation of (e, e'p) where the full spin dependence is included to a 
three dimensional numerical integration using the classical approximation. In like manner, 
we also replace the dependence in A by (rxp(r)')^. Finally, the approximate Coulomb 
distorted wavefunction of Eq. (19) with the A correction is given explicitly by [16] 

\]>(=‘=)(r) = ^Llrle±*'5((rxp'(r))2)gia(p'(r)-f)(rxp'(r))2gip'(r)-r^^^ ^20) 


with 5((rxp'(r))^) = 6o + & 2 (rxpX^"))^ + & 4 (rxp'(r))"^. Using this wavefunction and the hrst 
term of Eq. (5) we obtain the following four potential, which includes in an approximate 
way, the Coulomb distortion of the target nucleus 
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( 21 ) 


2l'‘(r) = 


4:71 


4pip /sin Y 


gi[5i((rxp'(r))2)+5j^((rxp'^(r))2)]^i(Ai-A/)gjq'(r)-r-^^/i^^ 


where A = a(rx p'(r))2(f ■p'{r)) and q'(r) = p'i{r) — p'j(r). The approximate potential 
of Eq. (21) is similar to the plane wave resnlt except for the phase factors and the radial 
dependence in the momentnm transfer. Unfortnnately, the spatial dependence in the phase 
factors makes a mnltipole decomposition of this potential impractical. However, since it 
is an analytical function it is straightforward to calculate interaction matrix elements by 
performing the three dimensional integration over r numerically. 


III. APPLICATION TO THE INCLUSIVE PROCESS 

For the inclusive cross section (e, e'), the longitudinal and transverse structure functions 
in Eq. (1) are bi-linear products of the Fourier transform of the components of the nuclear 
transition current density integrated over outgoing nucleon angles. Furthermore, it is the 
Dirac structure of the Moller potential which leads to the cross section containing one term 
with only longitudinal components of the current and a second containing only transverse 
components. However, the Dirac structure of the approximate potential in Eq. (21) is the 
same as the plane wave result. Therefore, even with Coulomb distribution included, albeit 
in an approximate way, the separation of the cross section into a longitudinal term and a 
transverse term persists. Explicitly, the structure functions for knocking out nucleons from 
a shell with angular momentum jh are given by 

St(9.u>)= E 0 ,0An f+ INyf'jdnp (23) 

^[^jb + i) J 

where the nucleon density of states pp = the 2 ;-axis is taken to be along q, and pb 

and sp are the z-components of the angular momentum of the bound and continuum state 
particles. The Fourier transfer of the nuclear current J^(r) is simply, 

N>^ = J (24) 
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The continuity equation has been used to eliminate the 2 ;-component (iV^) via the equation 
Nz = Note that when we use the approximate electron four potential along with 

current conservation to eliminate the z-component of the current we run into a problem since 
the momentum transfer q' depends on r both in magnitude and direction. In addition, the 
phase factors depend on r. To avoid generating additional terms we assume the direction 
of q'(r) is along the asymptotic momentum transfer q which defines the i-axis, and neglect 
the dependence on r in the phases and in q'(r), when taking the divergence of N. With 
this further approximation, current conservation implies uNq + q'(r)-N = 0. Using the 
approximate potential of Eq. (21), the cross section with for the inclusive reaction (e, e') can 
be written as 


^ w) + |tau^ I - ^]Sr(q ',»)} 


(25) 


and the transform of the transition nuclear current elements which appears in Sl and St 
are given by 


No = 



_^)2g%([rxp'(r)]2)gi5j-([rxp'^(r)]2)gi(Ai-A/)gjq'(r)-r 

g'(r) ^ 


Nt = 


j ( [r Xp' (r)] 2 )^*5/ ( [r Xp'^ (r)] 2 ) ^i( Ai -(r) -r 


JT{r)d^r 


(26) 

(27) 


As noted earlier, a multipole expansion of the approximate potential is not practical, and 
since the inclusive reaction (e, e') requires a sum over all occupied neutron and proton shells, 
numerical integration is very time consuming. In order to have a more practical procedure 
we choose to make additional approximations. Firstly, we neglect all of the phase factors(5 
and A) in Eqs. (26) and (27) but retain the r-dependence in q'. This returns us to the 
approximation we call LEMA. For light to medium nuclei this approximation is in good 
agreement with the full DWBA result regarding the shape as a function of energy transfer, 
but has a small discrepancy in magnitude. The magnitude is corrected with an overall factor 
of (AM ')2 cross section. However, for heavier nuclei, we noticed that LEMA with the 

magnitude factor was a very good approximation to the DWBA results for large electron 
scattering angle where the transverse term dominates, but deviated signihcantly for forward 
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electron angles where the longitudinal term has a signihcant contribution. Thus it appears 
that the phase factors play a signihcant role in the longitudinal term for large Coulomb 
distortion. 

We also noticed that for forward electron angles the low energy side of the DWBA 
quasielastic peak looks similar to the plane wave result. Thus, it appears that in the longi¬ 
tudinal term the phase factors are partially cancelling the effect of the effective momentum 
q'(r) on the low lo side of the quasielastic peak. We examined a number of simple ad — hoc 
modihcations to the longitudinal term as described by LEMA and the magnitude factor 
in order to reproduce the DWBA forward angle results for (e, e') reactions on in the 

quasielastic region. Based on a number of trials we propose the following “Fourier” transform 
for the charge component of the current 

^LEMA' ^ f eW'W'rj^(r)dV (28) 

Pi J 

where is the magnitude enhancement and q"(r) = pf(r) — Pj(r), p”{r) =p — 

^JgV(r')dr' and the factor A, which depends on the energy transfer to, is given by 
A = {u/ujqY with uq = Clearly for small a;, q" approaches the asymptotic value, 

while for ui past the quasielastic peak which is located approximately at o^o in fhe cross 
section, the effective momentum differs signihcantly from the asymptotic value. We have 
tested this ad — hoc prescription for 0.3<a;/ci;o<2.0 for a range of energies and nuclei and 
hud excellent agreement with the DWBA result. The transverse “Fourier” transform only 
contains the normalization factor and is given by 

j^lema' ^ f gWM-r (29) 

Pi J 

We will refer to the cross section for the inclusive reaction calculated with these two structure 
functions as the LEMA' result. Clearly and represent a modihed Fourier 

transform of the nuclear transition current. The approximation known as the EMA replaces 
q'(r) with q'(0) wherever it appears in Eqs. (26) and (27) for Nq and Nt and the phases 
are neglected as usual. We hud that for light nuclei the EMA is adequate, but it leads to 
large errors for nuclei as heavy as 
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In our analyses of quasielastic scattering, we use relativistic bound and continuum single 
particle wavefunctions. For the inclusive reaction we use continuum solutions for the out¬ 
going, but unobserved, nucleon which are in the same Hartree potential as the bound state 
orbitals [5]. This choice insures charge conservation and gauge invariance. Thus, in this 
relativistic single particle model the nuclear current matrix element is 

J^(r) = (30) 


where we use the free nucleon current operator 


> = + (31) 

2 m AT 

where /ip is the nucleon anomalous magnetic moment (for proton /ip = 1.793 and for neutron 
/ip = —1.91). The form factor Fi and F 2 are related to the electric and magnetic form factors 
Ge and Gm by 


Ge = 
Gm = 


Fi + 


4M2 ^ 


Fi -|- /ipF2 


(32) 

(33) 


We choose the standard result: 


Ge — Gm/{ft + — qI,/0-71) ^ (34) 

where in this formula qj, is in units of GeV^. 

In Fig. (4), we compare various approximations to the DWBA result as a function of the 
energy transfer u for two cases, incident electron energy Ei = 310 MeV, scattering angle 
6 *e = 143° and = 485 MeV, scattering angle 6^ = 60° data sets. The dotted line is the 
PWBA result, the diamonds are the DWBA result, the dash-dotted line is the EMA result, 
and the solid line is the LEMA' result. We notice that the EMA result is always lower than 
the DWBA result although the peak position is approximately in the right place. This lack 
of agreement with the EMA is not too surprising based on our previous examination of the 
wavefunction. Replacing the average value of the Coulomb potential between the origin and 
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the position r by the value at the origin is too large an error for r near the nuclear surface 
where most of the interaction takes place. Since the EMA is such a bad approximation to 
the full DWBA result, it should not be used as a basis for including the Coulomb phase 
terms in (e, e') or (e, e'p) reactions. Over the whole region, the LEMA' result is in excellent 
agreement with the DWBA result apart from the extreme wings of the quasielastic peak. 
The difference around the peak is less than 2 % and side parts are about 5 %. 



FIG. 4. The differential cross section for e') at two different electron energies and 

scattering angles. The dotted line is the PWBA result, the diamonds are the DWBA result, the 
dash-dotted line is the EMA result, and the solid line is the LEMA^ result. 

We illustrate our approximations for other kinematics by calculating the cross section 
at hxed momentum transfer q = 425 MeV/c with three different electron scattering angles 
{9e = 60°, 90°, and 143°) as shown Fig. (5). The EMA result is always lower than the DWBA 
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and the LEMA' results and is shifted toward large energy transfer by about 10 MeV. LEMA' 
again reproduces the DWBA cross sections quite well. 



IV. SEPARATION PROCEDURE AND STRUCTURE FUNCTIONS 

Since the LEMA' cross section for (e, e') has the same structure as the plane wave result, 
a Rosenbluth-type separation can be used to extract a longitudinal and a transverse contri¬ 
bution even in the presence of large Coulomb distortions. In fact, experimental evidence for 
this was given in one of the early papers on quasielastic scattering [4] where cross section 
measurements on ^38^ different electron scattering angles, but with the same en¬ 

ergy and momentum transfer, fell on a straight line when a Rosenbluth plot was made. Our 
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LEMA' approximation gives a theoretical explanation for this observation. Of course the 
separated structure functions are no longer bi-linear products of the simple Fourier trans¬ 
forms of the current components integrated over outgoing nucleon directions. Inclusion of 
Coulomb distortion within LEMA' results in an r-dependent Fourier momentum variable 
which differs for the longitudinal and transverse case. We have only been able check LEMA' 
for our particular model of the quasielastic process, so we can not prove that it applies to 
other models. However, based on previous work we know our model describes (e, e'p) and 
(e, e') quite well in the quasielastic region and its treatment of the spatial dependence in 
the nuclear charge and current is very realistic. Thus we believe that the Coulomb correc¬ 
tions that we calculate will be appropriate for any realistic nuclear model. We conclude 
that LEMA' is a good approximation for the inclusive cross section (e, e') in the quasielastic 
region. 

To illustrate the “Rosenbluth” separation we write S = Sl + xSt where S is the total 
structure function given by 

(35) 

Qfj. O'M hi ^edx! 

and X = [tan^^ - ^]/^. 

In Fig. (6), we compare the structure functions extracted using our calculations of (e, e') 
on 208 p^ using various treatments of Coulomb distortion for a momentum transfer of q = 425 
MeV/c and three different values of energy transfer around the peak of the cross section. In 
each case we hnd a very good ht to a straight line. Furthermore, we see that the intercept 
and slope extracted using the full DWBA and LEMA' calculations agree within 2% in all 
cases. The EMA and plane wave (PWBA) results clearly are in disagreement with the 
DWBA results. 
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FIG. 6. Rosenbluth separation plot of the cross section for ai q = 425 MeV/c. The 


solid line is the DWBA result, dashed line for the LEMA', the dotted line for the PWBA, and the 
dash-dotted line for the EMA. The unit of the total structure function is MeV“^. 



Given this result there are several ways to analyse experimental (e, e') data. If one 
has a model for the process under investigation, one can calculate the r-dependent Fourier 
transforms of Eqs. (28) and (29) and compare to the experimentally measured cross section. 
Or, one could make a Rosenbluth separation to obtain the LEMA' structure functions as 
a function of energy transfer uj and asymptotic momentum transfer q. Although these 
structure functions clearly have some dependence on electron energy and scattering angle 
in addition to their dependence on q and u, in our model of the quasielastic process this 
dependence is not very strong and we recommend that in comparing a theoretical model 
with the extracted LEMA' structure functions that electron energies and angles in the same 
range as used in the experiment be used. 

Lacking a model for the process being measured, one could take a simple model with 
suitable geometry and calculate the ratio of the cross section calculated with PWBA to one 
calculated with LEMA' and renormalize the measured cross section data by multiplying by 
this ratio. The resulting pseudo PWBA “data” could then be separated by a Rosenbluth 
procedure and would produce structure functions in terms of the Fourier transforms of the 
current components. Clearly this procedure introduces some error if the model amount of 
longitudinal and transverse contributions are not close to the amounts in the process being 
measured since LEMA' treats them somewhat differently. However, this difference is not 
too large and seems to be the best one can do. 

Finally one could assume that the theoretical model has the correct kinematics and 
spatial dependence, but that the magnitude of the longitudinal and/or transverse portions 
of the model is not correct. That is, one could use the model to calculate the longitudinal 
and transverse contributions to the cross section and then multiply each by a scale factor to 
be determined by making a least squares fit to the cross section data. 
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V. COMPARISON WITH EXPERIMENTAL DATA 


In Fig (7), we compare our theoretical results based on the relativistic “single particle” 
model [5] calculated with LEMA' to the Saclay data [8] for several electron angles {9f> = 
143°, 90°, 60° and 35°). The solid line is our model result for the cross section while the 
dotted line shows the longitudinal contribution to the cross section. Clearly the longitudinal 
contribution is quite small except for the forward electron scattering angle of 35° where it 
represents slightly over 50% of the cross section. Pion production is not included in our 
model, so the behavior at large oj is not expected to agree with the data. Clearly the 
agreement between the data and the calculation is not very good, which can be contrasted 
with quite good agreement between our calculations and the (e, e') data from Bates on 
"^°Ca [6]. Further, these may be a suggestion of longitudinal suppression if coupled with a 
transverse enhancement. 
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FIG. 7. The differential cross section for The solid line is the LEMA' result and 

the dotted line shows the longitudinal contribution. The data are from Ref. [8] 

We selected all of the Saclay data with energy transfer between 100 MeV and 200 MeV 
which had a clearly defined quasielastic peak in the cross section(218 data points), and used 
our model with LEMA' to calculate the longitudinal and transverse contributions to the cross 
section for each kinematical point in the data set. We performed a linear least squares fit 
by scale factors in front of the longitudinal and transverse contributions. This fit produced 
a factor in front of the longitudinal term of 0.69 and in front of the transverse term of 1.25. 
However, the fit is not very good since the data point is 60. The only conclusion 

we can make is that we do not find a 50% suppression of the longitudinal contribution in 
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“iospb^ Furthermore, it appears to us that the experimental data do not scale quite correctly. 
In Fig. (8) we show four sets of experimental data along with our model. Two of the sets 
are at forward angles with quite similar kinematics and two are at backward angles with 
similar kinematics. One of the forward angle data sets falls signihcantly below our model 
while the other is above. It would be difficult to modify the longitudinal and transverse 
strength in our model to get such a dramatic shift in behavior. For the larger angle case, the 
contribution of the longitudinal terms are negligible and again the kinematics do not vary 
so mnch between the two cases. The 310 MeV resnlt lies far above our calculation while the 
262 MeV result is only slightly above our calculation. 



u (MeV) « (MeV) 


FIG. 8. The inconsistent comparison with LEMA' and Saclay data. 
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VI. CONCLUSION 


We have developed a simple way of including Coulomb distortion in the Moller potential 
for inelastic electron scattering from medium and heavy nuclei. In this paper we have applied 
a simplihed version of this approximation to (e, e') reactions from medium and heavy nuclei 
in the quasielastic region. The previously used Effective Momentum Approximation (EMA) 
disagrees with the distorted wave analysis for nuclei as heavy as while the Local 

Effective Momentum Approximation with an enhancement factor and an ad —hoc correction 
to the longitudinal term (LEMA') reproduces the full DWBA calculation very well. Our most 
important Ending is that LEMA' allows the cross section to be separated into a longitudinal 
and a transverse contribution. However, the resulting structure functions depend on r- 
dependent Fourier transforms of the transition current components. We recommend several 
difierent procedures for using LEMA' for the analysis of experimental data including one 
where distortion eEects can be applied to the experimental data and “plane-wave” structure 
functions can then be extracted from the corrected data. 

We analysed the quasielastic (e, e') data from Saclay on using a relativistic single 

particle model and LEMA'. We do not find agreement with the data and even if we vary the 
longitudinal and transverse contributions in our model do not agree with the data. This is 
to be contrasted with our excellent agreement [6] with quasielastic data on and rough 
agreement with quasielastic data [4] on . We strongly recommend that additional (e, e') 
experiments on medium and heavy nuclei in the quasielastic region be carried out. The 
treatment of Coulomb corrections is no longer a serious hindrance to the analysis of such 
experiments. 
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